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ABSTRACT 

We present CHANDRA X-ray data of the nearby Magellanic Irregular dwarf star- 
burst galaxy NGC 4449. Contributions to the X-ray emission come from discrete 
point sources and extended diffuse emission. The extended emission has a com- 
plex morphology with an extent of ~ 2.4 x 1.6 kpc down to a flux density of 
1.3 x 10~ 13 erg s _1 cm -2 arcmin -2 . The best spectral fit to this emission is ob- 
tained with an absorbed, two temperature model giving temperatures for the two gas 
components of 0.28 ± 0.01 keV and 0.86 ± 0.04 keV, a total mass content of ~ 1O 7 M 
compared with a galactic mass of several 1O 1O M0 and a total thermal energy content 
of ~ 2.5 x 10 55 erg, with an average energy injection rate for the galaxy of a few 
10 41 erg s~ x . Comparison of the morphology of the diffuse X-ray emission with that 
of the observed Ha emission shows similarities in the two emissions. An expanding 
super-bubble is suggested by the presence of diffuse X-ray emission within what ap- 
pears to be a cavity in the Ha emission. The kinematics of this bubble suggest an 
expansion velocity of ~ 220 km s _1 and a mass injection rate of ~ 0.14Mq yr -1 , 
but the presence of NGC 4449's huge HI halo (r ~ 40 kpc) may prevent the ejection, 
into the inter-galactic medium (IGM), of the metal-enriched material and energy it 
contains. 

The arcsecond-resolution of CHANDRA has detected 24 X-ray point sources down 
to a completeness level corresponding to a flux of ~ 2 x 10~ 14 erg s _1 cm~ 2 , within 
the optical extent of NGC 4449 and analysis of their spectra has shown them to be 
from at least 3 different classes of object. As well as the known SNR in this galaxy, 
it also harbours several X-ray binaries and super-soft sources. The point source X- 
ray luminosity function, for the higher luminosity sources, has a slope of ~ —0.51, 
comparable to those of other starburst galaxies. 

Key words: ISM: jets and outflows - galaxies: individual: NGC 4449 - galaxies: 
starburst - X-rays: galaxies. 



1 INTRODUCTION 

NGC 4449 is a nearby Magellanic Dwarf Irregular starburst 
galaxy which has an inclination of 56.2° and a mass of 
~ 4 x 10 10 Mq (Bajaja et al. 1994 - mass has been scaled 
for our assumed distance which is discussed later). As such 
it allows another opportunity for the study of the starburst 
phenomenon and its effect on galaxy evolution in the lo- 
cal Universe. Dwarf galaxies as the basic building blocks 
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in the hierarchical merging cosmology scenario are likely to 
have harboured the earliest sites of star-formation in the 
Universe and so their study in the local Universe can give 
insight into the evolution of such objects at high redshift. 
Observations of local edge-on starburst galaxies (Strickland 
et al. 2000; Weaver 2001) are presenting a picture of kpc- 
scale, soft X-ray emitting, bi-polar outflows in the form of 
galactic winds transporting mass, newly synthesized heavy 
elements and energy into the IGM. These winds result from 
the pressure driven outflows along these galaxy's minor axis 
produced from the efficient thermalization of the mechani- 
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cal energy from the supernovae (SN) explosions and stellar 
winds of the massive stars in their OB associations and su- 
per star clusters (SSC). Inclined galaxies such as NGC 4449 
containing such phenomena present a less clearly observable 
picture. However, observations of how absorption and tem- 
perature of the diffuse X-ray emission varies within them 
and links between Ha, HI and X-ray morphology suggest a 
similar scenario also applies to these objects. 

NGC 4449 has been observed across the electromagnetic 
spectrum and displays some both interesting and unusual 
characteristics. Radio observations have shown it to contain 
an extremely luminous SNR (Bignell & Seaquist 1983) , have 
a very extended HI halo (~ 40 kpc in radius) which appears 
to be rotating in the opposite direction to the gas in the core 
of the galaxy (Bajaja et al. 1994) and to have a large-scale 
ordered magnetic field (Klein et al., 1996). The angular ve- 
locity of the rotation of the gas in the centre of the galaxy is 
low (~ 18 km s" 1 Hunter et al. 1998), suggesting a low es- 
cape velocity but the huge HI halo may prevent escape of the 
hot ejecta of the starburst region from the galaxy as a whole. 
The galaxy has both a high star-formation rate (~ 0.2 Mq 
yr _1 , Thronson Jr. et al. 1987) and supernova rate, contains 
numerous star clusters (~ 60, Gelatt et al. 2001) with the 
central one appearing to be young (~ 6 — 10 Myr Boker et 
al. 2001) and a spherical distribution of older stars with a 
mean age of 3 — 5 Gyr (Bothun 1986). It has also been shown 
to contain molecular cloud complexes from CO observations 
(Hunter & Thronson Jr. 1996) and has an infra-red luminos- 
ity (10 - 150 /Lim) of 3.7 x 10 42 erg s" 1 (Thronson Jr. et al. 
1987, corrected for the distance assumed here - see below). 
The ionised gas it contains shows a very disturbed morphol- 
ogy and includes what appear to be many Ha bubbles, shells 
and filaments (Hunter & Gallagher 1990; 1997). The kine- 
matics of the HII regions within the galaxy are very chaotic 
and could be the aftermath of a collision or merger (Hart- 
mann et al. 1986). There is a gas-rich companion galaxy 
lying at a projected distance of ~ 40 kpc, coincident in pro- 
jection with the outer edge of the extended HI halo, which 
could have been involved. Previous X-ray observations of 
this galaxy have detected the presence of discrete emission 
from point sources and diffuse emission from hot gas (Delia 
Ceca et al. 1997; Vogler & Pietsch 1997). 

The distance to NGC 4449 is not well established and 
ranges from 2.93 Mpc (Karachentsev & Drozdovsky 1998) 
to 5 Mpc (Aaronson & Mould 1983). We have chosen to 
adopt the lower value throughout this analysis as optical 
observations are showing objects within NGC 4449 to be 
resolved (Whiting 2002, private communication) suggesting 
a lower distance estimate to be more appropriate. 

In Section 2 we describe the CHANDRA observation. 
The X-ray emission from the point sources is discussed in 
Section 3 and that from the diffuse emission is considered 
in Section 4. Section 5 contains a more general discussion 
of the relationship of the X-ray emission to that from other 
wavebands, the affects of the out-flow of the X-ray emission 
on NGC 4449 and the morphology of the X-ray emission, 
whilst our main conclusions are summarized in Section 6. 



2 OBSERVATIONS AND ANALYSIS 

A 30 ks observation of NGC 4449 was obtained on Feb 
4" 1 — 5 th 2001. Analysis was carried out on the data con- 
tained within the S3 chip of the ACIS— S instrument on 
board CHANDRA using CIAO (version 2.2.1), HEASOFT 
(version 5.1), XSPEC (version 11.1.0) and ASTERIX (ver- 
sion 2.3-bl). The data were reprocessed (CALDB version 
2.9) and filtered to remove periods of flaring and lower than 
average count rates using the 3<r clipping technique of the 
lc_clean.sl script (available from the CHANDRA website). 
After these processes were completed, a total of 22.12 ks 
of useful data remained. The data were then further fil- 
tered to contain data in the energy band from 0.3 — 8.0 keV 
(amounting to a total of ~ 22300 counts on the S3 chip, 
of which ~ 10700 come from within the D25 ellipse and 
of these ~ 6100 are due to diffuse emission) and were also 
background subtracted using the appropriately scaled CXC 
background event data set. (The one utilized for this obser- 
vation had a total exposure of 137391s.) After production of 
an exposure map, the point sources present in the data were 
detected using the wavdetect tool, contained in the CIAO 
software, run with an exposure map, default wavlet scales 
of 2.0 and 4.0 pixels and a source significance threshold of 
5 x 10 -7 , a value of ~ (number of pixels) -1 which should 
have limited the number of false detections to ~ 1. 32 point 
sources were detected and these are shown in Fig. ^overlaid 
on a smoothed image of the S3 chip (smoothed with a Gaus- 
sian having a FWHM of 4 pixels ~ 2"), while their positions 
and background subtracted count rates are listed in Table 1. 
Of these 32 sources, 8 do not lie within the optical extent of 
the galaxy, as measured by the D25 ellipse (de Vaucouleurs 
et al. 1991), also shown in Fig-E] The 24 sources within the 
optical extent of the galaxy are the ones most likely to be as- 
sociated with NGC 4449 and so are the only ones considered 
further in this analysis. In comparison, when the data were 
point source searched using celldetect with an exposure map 
and default settings, only 18 sources were detected on the 
S3 chip. It was obvious from visual inspection of the data 
that several significant sources were missed by celldetect. In 
addition on comparison with expected source counts from, 
'The CHANDRA Deep Field North Survey' (Brandt et al. 
2001), we would expect ~ 16 — 20 background sources on 
the S3-chip at our level of completeness (see section 3.1 for 
further details). For these reasons, we would expect more 
than the 18 sources detected by celldetect and so the results 
from wavdetect have been adopted as being more reliable. 

After detection, the point sources were subtracted from 
the data so that the diffuse X-ray emission contained within 
NGC 4449 could be investigated. The sources were simply 
blanked, using areas equivalent to the 3a detection ellipses 
from wavdetect and allowance was made for the lost flux 
from these areas, which amounted to only 0.5% of the total 
diffuse emission, in later analysis. Also, a 3 colour adap- 
tively smoothed image of the galaxy has been produced us- 
ing csmooth in fft (fast fourier transform) mode and this 
is shown in Fig. |5] The lower and upper sigmas for deriv- 
ing the smoothing kernel were set at 2 and 5 respectively. 
This image shows that some of the sources are very soft and 
that the diffuse emission appears to consist of components of 
varying temperatures. The red, green and blue images used 
correspond to energy bands of 0.3 — 0.8 keV, 0.8 — 2.0 keV 
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Figure 1. Low resolution (smoothed using a Gaussian with FWHM of 4 pixels ~ 2 ), background subtracted image of the S3 chip 
field of view marked with the 32 point sources listed in Table 1 and the D25 ellipse, representing the optical extent of NGC 4449. The 
extended diffuse emission associated with the galaxy is clearly visible. 



and 2.0 — 8.0 keV respectively. These bands were chosen so 
that the total number of counts in each band was approx- 
imately equal thereby maximizing the S/N ratio for the 3 
bands and allowing point source searching of these individ- 
ual bands to be performed most effectively. These bands 
are also in line with those used in the analysis of other 
starburst galaxies (e.g. NGC 253, Strickland et al. 2002). 
The diffuse emission is seen to be extended and the size of 
the region it occupies extends 2.4 kpc from NNE to SSW 
and 1.6 kpc from WNW to ESE down to a flux density of 
1.3 x 10 -13 erg s" 1 cm" 2 arcmin -2 . 

The background subtracted spectrum of the total X- 
ray emission from within the D25 ellipse of NGC 4449 was 
fitted using the modified Levenberg-Marquardt method and 
standard \ 2 statistic from XSPEC, with an absorbed 2 ther- 
mal component plus power law fit. The resulting spectral 
fit is shown in Fig. [3] The absorbing column density was 
considered to be two components, one due to the inter- 
stellar medium (ISM) within the Milky Way (Galactic col- 
umn density, Nh = 1-4 x 10 20 cm" 2 ) and the other due 
to the ISM of NGC 4449. The 2 absorbed thermal com- 
ponents were modelled using the wabs (photo-electric ab- 
sorption using Wisconsin cross-sections) and mekal thermal 



plasma codes within XSPEC, initially Galactic absorption 
and an abundance of 0.276 Solar (the value obtained by 
Martin, 1997, for NGC 4449) were assumed. Two thermal 
plasmas were used to represent the diffuse fraction of the 
total X-ray emission as images in the soft (0.3 — 0.8 keV) 
and medium(0.8 — 2.0 keV) energy bands, see Fig. [5] show 
different spatial distributions of the hot gas and also, multi- 
phase models are seen to be needed to best fit the emission 
from other starburst galaxies (e.g. NGC 253, Strickland et 
al. 2002; NGC 1569, Martin et al. 2002). A single tempera- 
ture fit to NGC 4449 also gave a less robust statistical fit and 
in particular when fitted to just the diffuse emission does not 
have a wide enough energy distribution to fit the spread seen 
in the data. The two-temperature model gives an absorption 
corrected flux of (2.39±g;i|) x 10~ 12 erg s" 1 cm" 2 , which 
corresponds to a total X-ray luminosity, in the 0.3 — 8.0 keV 
band, of (2AQ±lf Q ) x 10 39 erg s" 1 . The fitted column density 
for NGC 4449 was (1.31 ± 0.52) x 10 21 cm" 2 , with the two 
thermal components having temperatures of 0.27 ±0.01 keV 
and 1.01T0.06 keV and a fitted abundance of 0.29±g;g| Z . 
The power law component had a fitted photon index of 
F = 2.05 ± 0.46. Of this total emission, ~ 60% is from the 
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Figure 2. Adaptively smoothed 3 colour image of NGC 4449 (red: 0.3 - 0.8 keV, green: 0.8 - 2.0 keV and blue: 2.0 - 8.0 keV). The 
extent of the diffuse emission is ~ 2.4 kpc from NNE to SSW and ~ 1.6 kpc from WNW to ESE. The variation of hardness in the spectra 
of the resolved sources is evident from their varying colours and is indicative of the presence of different types of sources, such as XRBs, 
SNRs and SSS. 



resolved point sources, ~ 30% is from the cooler thermal 
component and ~ 10% is from the hotter component. 



3 POINT SOURCES 
3.1 Source Spectra 

A background subtracted spectrum was extracted for each 
of the 24 sources within the D25 ellipse. (The background 
spectra were taken from the aforementioned background file 
and were appropriately scaled. The use of the background 
files meant that the background spectrum was taken from 
the identical position on the CCD and avoided the prob- 
lem of contamination from adjacent sources.) The spectra 
were grouped so that each had a minimum of 5 data bins 
and when possible a minimum of 10 counts per bin. No at- 
tempt was made initially to fit the data where a source had 
less than 50 counts after background subtraction. The fits 
obtained for the 5 sources with the highest count rates are 
shown in Fig^] Comparison of the low count spectra to those 
fitted allowed a rough fit to be made to these sources indi- 
vidually by using the models of sources with similar spectra 
and hardness ratios (see below). The best fits obtained for 
the sources are summarized in Table 2. The same start pa- 
rameters were assumed for column densities and abundances 
as were used for fitting the total emission. 

For each of the individual source fits, the unabsorbed 
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channel energy (keV) 



Figure 3. Fitted spectrum of the total X-ray emission from 
NGC 4449. The fit shown is an absorbed two temperature plus 
power law fit, with a column density of Njj = 1.40 X 10 20 cm" 2 
assumed for Galactic absorption and a fitted column density 
of N H = (1.31 ± 0.52) X 10 21 cm" 2 for the absorption within 
NGC 4449. The fitted temperatures of the two thermal compo- 
nents are 0.27 ± 0.01 and 1.01 ± 0.06 keV with the fitted abun- 
dance for both components being (0.29 ±o'g| Zq and the fit- 
ted photon index of the power law component having a value 
of r = 2.05 ± 0.46. The individual contributions of the separate 
components of the fit are also included. 
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Table 1. Positions and count rates of the 32 sources detected in 
the NGC 4449 CHANDRA S3 chip data. Column 1 is the source 
number ordered in increasing R.A.. Columns 2 and 3 give the 
R.A. and Dec. of each source and Column 4 lists their background 
subtracted count rates. 



Source RA (h m s) Dec (° ' ") Count Rate 

(xlCT 3 cts s" 1 ) 



1 


12 


27 


44.24 


44 


07 


23. 


.0 


0.13 ±0.08 


2 


12 


28 


00.65 


44 


04 


32 


.5 


0.37 ±0.14 


3 


12 


28 


01.36 


11 


05 


29. 


.9 


0.63 ±0.18 


4 


12 


28 


02.07 


11 


08 


08. 


.0 


5.09 ±0.48 


5 


12 


28 


03.89 


14 


05 


43. 


.8 


5.36 ±0.50 


6 


12 


28 


06.82 


41 


05 


28. 


.4 


0.83 ±0.21 


7 


12 


28 


07.13 


41 


04 


14. 


.4 


0.96 ±0.22 


8 


12 


28 


07.29 


41 


07 


21. 


.6 


2.96 ±0.37 


9 


12 


28 


07.34 


44 


04 


53. 


.8 


5.65 ±0.52 


10 


12 


28 


09.32 


44 


05 


08. 


.3 


43.18 ± 1.42 


11 


12 


28 


09.71 


44 


05 


53. 


.1 


7.53 ±0.59 


12 


12 


28 


10.40 


41 


05 


58. 


.1 


1.96 ±0.32 


13 


12 


28 


10.52 


11 


01 


59 


.5 


1.04 ±0.24 


1 1 


1 9 


9vi 


iu.yo 


AA 
11 


n 

uo 


■ > O . 


.u 




15 


12 


28 


10.95 


44 


06 


48. 


.5 


38.60 ± 1.33 


16 


12 


28 


10.97 


44 


04 


59 


.5 


0.94 ±0.23 


17 


12 


28 


11.20 


11 


06 


38 


.5 


4.13 ±0.44 


18 


12 


28 


11.33 


41 


05 


37 


.0 


0.50 ±0.19 


19 


12 


28 


11.84 


11 


05 


30 


.0 


1.80 ±0.32 


20 


12 


28 


11.98 


11 


06 


41 


.5 


6.50 ±0.56 


21 


12 


28 


12.02 


11 


05 


58. 


.4 


2.18 ±0.33 


22 


12 


28 


13.00 


44 


08 


44. 


.1 


0.35 ±0.13 


23 


12 


28 


13.26 


41 


06 


55. 


.8 


9.61 ±0.68 


24 


12 


28 


14.93 


44 


04 


31. 


.8 


0.36 ±0.14 


25 


12 


28 


15.64 


41 


05 


37 


.1 


0.63 ±0.19 


26 


12 


28 


16.41 


41 


07 


39. 


.1 


0.42 ±0.14 


27 


12 


28 


17.83 


41 


06 


33. 


.9 


51.35 ± 1.55 


28 


12 


28 


19.02 


44 


05 


44. 


.6 


1.69 ±0.28 


29 


12 


28 


22.44 


44 


03 


39. 


.7 


0.39 ±0.14 


30 


12 


28 


25.60 


41 


04 


37. 


.9 


1.41 ±0.26 


31 


12 


28 


31.86 


11 


02 


42. 


.5 


1.74 ±0.33 


32 


12 


28 


35.02 


11 


02 


56. 


.0 


0.26 ±0.12 



flux in the 0.3 — 8.0 keV energy band and luminosity (assum- 
ing a distance of 2.93 Mpc for NGC 4449) were calculated 
and these luminosities are also shown in Table 2 along with 
the errors determined from the 90% confidence regions for 
each model's normalization. In addition the data was point 
source searched in the three different energy bands used to 
produce the three colour image, (soft 0.3 — 0.8 keV, medium 
0.8 - 2.0 keV and hard 2.0 - 8.0 keV) so that hardness ra- 
tios for the individual sources within the D25 ellipse could be 
calculated for the sources detected in more than one of the 
energy bands. These ratios are shown in Table 3. The cumu- 
lative Log(N)-Log(Lx) plot for the 24 fitted sources possibly 
associated with NGC 4449 is shown in Fig. This has been 
fitted with a single power law with a slope of ~ —0.51 for the 
higher luminosity sources. The highest luminosity source is 
not included in the fit. To do so results in an offset between 
the data and the fitted line but does not alter the slope of 
the fit. The spectral fit to this source (source 28) does sug- 
gest it is a very unusual source as it is very soft, lies behind 
a high column density and has an extremely high luminos- 
ity (> 10 41 erg s~ 4 ). The fact that its spectrum contains 
< 40 counts in total brings into question the robustness of 



the fit and for these reasons it has been excluded from the 
fit. Sources with luminosities below an absorption corrected 
luminosity of 2.18 x 10 37 erg s _1 , equivalent to a flux of 
2.12 x 10 -14 erg s _1 cm -2 , were also not included in the 
fit and these figures represent the completeness limit for the 
data. The fitted value for the slope is comparable to the 
same fits for M82 and the Antennae where the slope was 
determined to be - -0.45 ± 0.06 (Zezas et al. 2001). The 
absorption corrected contribution to the X-ray luminosity 
of NGC 4449 from unresolved point sources is estimated to 
be ~ (5.62±t;§§) x 10 37 erg s"\ equivalent to ~ 5% of the 
diffuse emission (see section 4). 

Earlier ROSAT HRI observations of this galaxy (Vogler 
& Pietsch 1997) detected only 7 point sources within the D25 
ellipse. Of these 7, at least 3 of them have been resolved into 
more than one source using the superior arcsecond spatial 
resolution of CHANDRA and an additional 13 sources have 
also been detected. The source identified as the supernova 
remnant (SNR), (X4) in the ROSAT data corresponds to 
source 15 in our CHANDRA data when the source position 
is compared to the accurate position for the SNR obtained 
by Bignell & Seaquist (1983), from radio measurements and 
there are also 2 additional sources resolved at the position of 
X4 in the ROSAT data. This highlights the fact that lumi- 
nosity estimates based on pre-CHANDRA data are effected 
by source confusion. Also of note is the fact that source 27 
in our data ( X7 in the ROSAT data) is much more lumi- 
nous than it was in 1994 and now has the highest count-rate 
for a single source in the galaxy rather than being fourth 
highest. See Table 4 for a comparison of the 0.1 — 2.4 keV 
absorption corrected luminosities of sources that were both 
detected by ROSAT and not resolved into more than one 
source by CHANDRA. 

In general, the sources seem to split into four classes 
which can be represented by sources 14, 15, 27 and 2. Sources 
like source 14 are very soft, luminous sources (super-soft 
sources, SSS) which show little emission above 1 keV. Source 
15 is the SNR and at least 5 other sources show similar spec- 
tra and/or hardness ratios. Sources like source 2 have spec- 
tra that contain hard emission up to ~ 5 keV and sources in 
this group seem to be best fitted by power law models. They 
are generally fitted with a higher column density and lower 
photon index than other sources and may be background 
AGN. Sources like source 27 also have harder components 
to their spectra, showing emission above 5 keV and are fit- 
ted by power laws with 1.3 < F < 3.0, typical of X-ray 
binaries (XRB). See Section 3.2 below for further discus- 
sion on the brightest sources. Fig. shows a comparison of 
the hard and soft hardness ratios for the sources detected 
in all 3 energy bands. The hard ratio is (h — m)/(h + m) 
while the soft ratio is (m — s)/ (m ± s), where s is the count 
rate in the soft band (0.3 — 0.8 keV), m is the count rate 
in the medium band (0.8 — 2.0 keV) and h is the count rate 
in the hard band (2.0 — 8.0 keV). The sources represented 
by squares were only detected in 2 of the 3 energy bands 
and so their positions represent the most extreme positions 
possible for them on this plot. Those at the bottom need to 
move vertically upwards while those to the right will need to 
move horizontally to the left. The sources appear to clump 
in groups on this plot with the SNR and SSS sources being 
softer than those we have classified as XRB and AGN. The 
right hand plot in this figure shows some theoretical tracks of 
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Table 2. Best fitting single component models for the spectra of the 24 point sources. Column 1 gives the source numbers as shown on 
Fig-E Column 2 contains the best fit absorbed single component models obtained. In each case, the absorbing column due to the Milky 
Way, wabs Ga ; is assumed to be Njj = 1.4 X 10 20 cm -2 Column 3 is the column density obtained from fitting the wabs (photo-electric 
absorption using Wisconsin cross-sections) component, column 4 the fitted temperature for the mekal, brems (thermal bremsstrahlung), 
bb (black body), and nei (non-equilibrium ionization) models for the thermal components, column 5 is the photon index for the po 
(power law) models, column 6 the statistic of each fit and column 7 the absorption corrected luminosity for each source. The errors 
shown on the luminosities are from the 90% confidence regions (1.64cr) for each model's normalization. Where only one value is shown 
this is an upper limit from the same 90% confidence regions. Column 8 lists the type of source each object might be from comparison 
with the spectra and hardness ratios of the identified sources within the galaxy. 



Source Model Njj (cm 2 ) kT (keV) Y x 2 /d.o.f Luminosity (erg s *) Source Type 



2 


wabs Gai (wabs(po)) 


3.14 


X 


10 21 




0.61 


1.78/5 


(7.10± 731 ) 


X 


10 36 


AGN 


3 


wabs Ga ; (wabs(mckal)) 


3.68 


X 


10 21 


0.90 




0.49/2 


(5.87dz|;|f ) 

l c 1 c_i_0 75 A 


X 


10 36 


SNR/XRB 


5 


wabs Ga( (wabs(po)) 


6.58 


X 






1.30 


1.39/8 


X 


10 37 


XRB 


6 


wabs Ga ( (wabs(mckal)) 


1.31 


X 


10 21 


0.82 




2.01/2 


(5.51± 3 - 52 ) 


X 


10 36 


SNR/XRB 


7 


wabs Ga! (wabs(po)) 


1.39 


X 


10 21 




0.40 


1.00/2 


(2.18± - 98 ) 


X 


10 37 


AGN 


8 


wabs Ga! (wabs(po)) 


1.53 


X 


10 20 




1.33 


0.38/3 




X 


10 37 


XRB 


9 


wabs Gai (wabs(po)) 


3.93 


X 


10 20 




2.88 


7.87/8 




X 


10 37 


SNR/XRB 


10 


wabs Gai (wabs(brcms)) 


9.29 


X 


10 20 


1.42 




83.37/69 


(2 24+°' 12 i 


X 


10 38 


SNR/XRB 


11 


wabs Gai (wabs(po)) 


3.50 


X 


10 21 




1.64 


9.07/13 


(8.91±i; 2 3) 


X 


10 37 


XRB 


12 


wabs G(l ; (wabs(bb)) 


4.15 


X 


10 20 


0.091 




4.46/4 


(1.28± 2 - n ) 


X 


10 37 


SSS 


11 


wabs Ga ; (wabs(bb)) 


8.96 


X 
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X 


10 3S 


sss 


15 
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1.26 


X 
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2.58 




65.50/60 


(o Q4+0.14\ 


X 


10 38 


SNR 


16 


wabs Gai (wabs(po)) 


1.40 


X 


10 22 




2.53 


0.90/3 


(3.i4±;-i|) 


X 


10 37 


XRB 


17 


wabs Ga ; (wabs(nei)) 


7.45 


X 


10 20 


19.21 
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C2 c;i+o:50' l 


X 


10 37 


SNR/XRB 


IS 


wabs Ga j (wabs(nei)) 


3.82 


X 


10 20 


3.05 




1.03/3 


(3.02+. 1 - 37 ) 


X 


10 36 


SNR/XRB 


19 


wabs Gai (wabs(po)) 


1.31 


X 


10 21 




1.35 


0.47/3 


(o S3+ 1 ' 55 ") 


X 


10 37 


XRB 


20 


wabs Ga! (wabs(po)) 


3.45 


X 


10 21 




2.76 


8.20/10 


(7.43±i;0|) 


X 


10 37 


SNR 


21 


wabs Ga ; (wabs(bb)) 


2.45 


X 


10 21 


0.080 




1.71/2 


(7.19±i-°7) 


X 


10 37 


SNR/SSS 


23 


wabs Ga ; (wabs(nei)) 


9.00 


X 


10 21 


1.30 




12.07/16 


(8.06±;;*f) 


X 


10 38 


SNR/XRB 


24 


wabs Ga ; (wabs(po)) 


1.10 


X 


10 22 




0.97 


2.25/6 


(5.79± 5 - 22 ) 


X 


10 36 


AGN 


25 


wabs G(1 ; (wabs(bb)) 


1.11 


X 


10 21 


0.130 




0.62/2 


(5.17± 107 - 71 


) x 10 36 


SSS/XRB 


26 


wabs Ga j (wabs(brcms)) 


4.72 


X 


10 21 


54.85 




0.74/2 


(4.99± 357 ) 


x 10 36 


AGN 


27 


wabs Ga ; (wabs(po)) 


6.49 


X 


10 21 




1.91 


64.38/86 


(a iq+0.42\ 


X 


10 38 


XRB 


28 


wabs G(l ; (wabs(bb)) 


5.01 


X 


10 21 


0.050 




1.28/3 




x 10 41 


SSS 



how a sources position would vary depending on which type 
of spectral model is assumed for it, what value the absorbing 
column density that it lies behind has and what tempera- 
ture or photon index is assumed for the model, for compari- 
son with the plotted source positions. An attempt has been 
made to classify the sources using the shape and range of 
spectra and their positions on the hardness ratio plot. The 
suggested classifications are shown in column 8 of Table 2. 
The D25 ellipse covers ~ 25% of the S3-chip and as such 
it would be expected from the background number counts 
quoted in section 2 that ~ 4 — 5 background sources should 
be found within this area. This figure agrees well with the 
number of sources classified as potential background AGN. 
In contrast to this, the 8 point sources detected outside the 
D25 ellipse are only about a half to two-thirds of what would 
be expected, suggesting that our completeness level estimate 
may be at slightly too low a flux. 

3.2 More Detail on the Brightest Sources 

3.2.1 Source 10 - SNR? 

This source lies in the soft region of the hardness plot of 
Fig. ©to the left of the area occupied by the known SNR. It 
is best fitted by an absorbed thermal Bremsstrahlung model. 
Its position would also suggest it should have a low column 
density as it lies to the left of the theoretical tracks. This 



does not appear to be the case and fitting with other models 
gives an even higher column density for a power law model or 
unphysical results. This source was also detected (source X7) 
in the earlier ROSAT observations and shows no evidence 
for variability as seen in Table 4. Hence, it seems more likely 
that it is a SNR than an XRB. 



3.2.2 Source U - SSS 

This source is very soft and it can be fitted in a similar way 
to CAL 87 (Ebisawa et al. 2001) using an absorbed black- 
body component and several absorption edges, indicative of 
a white dwarf binary system. For such a fit, three absorp- 
tion edges are included at energies of 0.40 keV, 0.59 keV 
and 0.85 keV the value of Nh increases to 3.24 x 10 21 cm -2 
and kT decreases to 0.065 keV, compared to the one com- 
ponent fit shown in Table 2, while the absorption cor- 
rected luminosity in the 0.3 — 8.0 keV band increases to 
~ 4.05 x 10 39 erg s _1 . The absorption edges are possi- 
bly connected with the presence of the following elements: 
0.40 keV, NVI or CVI; 0.59 keV, OVII or NVII; 0.85 keV, 
OVIII and/or OVII (the low source counts making a defi- 
nite identification difficult). The blackbody temperature fit- 
ted here lies in the range (T ~ 60 — 80 eV) where both 
OVII and OVIII edges would be expected to be present with 
the latter being more abundant (Ebisawa et al. 2001). This 
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Figure 4. Spectra and best fitting models for the brightest sources detected within the D25 ellipse. Results are shown for source 10 (top 
left), source 14 (top right), source 15 (middle left), source 23 (middle right) and source 27 (bottom). All the fits are for the absorbed one 
component models shown in Table 2. More detailed fits for the brighter sources are discussed in Section 3.2 of the text. 



source has a greatly increased absorption corrected luminos- 
ity compared with the earlier R OS AT observations however 
the difference is in part due to the different spectral models 
used to fit the spectrum. The improved spectral resolution 
of CHANDRA compared with ROSAT clearly shows that 
a 5 keV thermal bremsstrahlung model is inappropriate for 
this source, as all its emission lies below 1 keV. In Table 4, 
the absorption corrected luminosity is also much higher than 
that shown in Table 2 and quoted above due to the value 
in Table 4 including emission down to 0.1 keV rather than 
0.3 keV. These facts make it difficult to assess if there has 
been any intrinsic increase in the sources luminosity over 
the last decade. When the CHANDRA data are fitted with 



the 5 keV thermal Bremsstrahlung model there is agreement 
with the ROSAT HRI data, within errors. 

3.2.3 Source 15 - SNR 

This source has been identified as a young and very luminous 
SNR embedded in an HII region from observations in sev- 
eral different wavebands, radio (Seaquist & Bignell 1978), 
optical (Balick & Heckman 1978, Blair et al. 1983), UV 
(Blair et al. 1984) and X-ray (Blair et al. 1983; Vogler & 
Pietsch 1997). Fitting this source with an absorbed, non- 
equilibrium ionization model (as used for modelling other 
SNR, e.g. Yokogawa et al., 2002) gives an X-ray temper- 
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Table 3. Hardness ratios, where possible, for the 24 sources detected within the D25 ellipse. Soft band, s, 0.3 — 0.8 keV, medium band, 
m, 0.8 — 2.0 keV and hard band, h, 2.0 — 8.0 keV. Where no counts are shown, wavdetect failed to detect the object in that energy band. 
Column 1 gives the source numbers as shown on Fig. Columns 2 — 4 are the counts in the 3 different energy bands and columns 5 and 
6 give the values of the hardness ratios calculated as detailed in the column headings. 



Source 


Counts in 


Counts in 


Counts in 


(m - s) 


(h - m) 




Soft Band 


Medium Band 


Hard Band 


(m + s) 


(h + m) 


2 




2.0 ± 1.4 


3.9 ± 2.0 




n 32 + 43 


3 




8.9 ± 8.0 








5 


25.7 ±5.1 


58.7 ± 7.7 


34.7 ± 5.9 


n 30 + n 1 2 


_n 26 + 1 1 

\J.£i\J _1_ U.J-J- 


5 




8.9 ± 8.0 








7 




3.9 ± 2.0 


7.9 ± 2.8 




n 34 _i_ n 31 


8 


1 2 8 + 3 fi 


31.8 ± 5.7 


20 s + 4 fi 

.iiL/.O _1_ rfc.U 


n 43 + 1 7 

W.rtO _1_ W.J- 1 


—0 21 +014 


9 


w 2 + 7 5 


57.4 ± 7.6 


6.9 ± 2.6 


n 02 + 1 


_n 7Q + n Ifi 


10 


KQQ 7 _|_ 23 2 


Q'JK nxisi 

JOU.i? _1_ J-O.Tl 


03 6 + 100 


_n 23 + 03 


_n ^7 + n 06 

U.iJl _1_ u.uu 


11 


23.9 ± 5.0 


88.0 ± 9.4 


54.6 ± 7.4 


0.57 ±0.11 


-0.23 ±0.09 


12 


39.6 ±6.4 










11 


349.0 ± 18.7 


7.9 ± 2.8 




-0.96 ±0.07 




15 


141.7 ± 12.0 


590.5 ±24.4 


126.5 ± 11.3 


0.61 ±0.04 


-0.64 ±0.11 


16 




10.5 ±3.3 


6.9 ± 2.6 




-0.21 ±0.17 


17 


30.1 ± 5.6 


54.1 ± 7.4 




0.29 ±0.49 




18 












19 




22.4 ±4.9 


14.7 ±3.9 




-0.21 ±0.17 


20 


30.1 ± 5.6 


91.1 ±9.6 


19.7 ±4.5 


0.50 ±0.10 


-0.64 ±0.11 


21 


14.8 ±4.0 


31.7 ± 5.7 




0.36 ±0.16 




23 


51.8 ± 7.3 


129.6 ± 11.4 


27.6 ±5.3 


0.43 ± 0.09 


-0.65 ±0.09 


24 






4.9 ±2.2 






25 












26 






4.9 ±2.2 






27 


32.4 ±5.7 


669.5 ±25.9 


448.7 ±21.2 


0.91 ±0.05 


-0.20 ±0.03 


28 


36.6 ±6.1 











Table 4. Comparison of the absorption corrected luminosities for sources 10, 14, 23 and 27 obtained from ROSAT PSPC observations 
(11/1991), ROSAT HPJ observations (12/1994) and CHANDRA ACIS-S observations (02/2001). The ROSAT values are from Vogler 
& Pietsch (1997) and are for 5 keV thermal bremsstrahlung spectra, corrected for Galactic absorption. The first column of CHANDRA 
values are for the best fits discussed in the text, corrected for Galactic absorption. All luminosities are for the 0.1 — 2.4 keV energy band 
and the CHANDRA values have been scaled to account for the differences in assumed distances to NGC 4449. The errors shown on the 
CHANDRA data are from the 90% confidence levels (1.64cr) of the normalizations of the fits. Columns 1 and 2 are the source numbers 
allocated in the respective data sets and columns 3 — 5 give the absorption corrected luminosities fitted to the 3 data sets. Column 6 
shows the absorption corrected luminosities obtained for the CHANDRA data when they are fitted with a 5 keV thermal bremsstrahlung 
model, corrected for Galactic absorption. 



Chandra 


ROSAT 


Lx(PSPC) 
(xlO 38 erg s" 1 ) 


L X (HRI) 
(xlO 38 erg s _1 ) 


ix(ACIS-S) 
(xlO 38 erg s -1 ) 


L x (ACIS-S 5 keV) 
(xlO 38 erg s _1 ) 


10 


XI 


2.92 ±0.43 


3.10 ±0.37 


3 24+0- 18 
105.45±? 53 8 
14.60±i;£ 


2.58±»;« 


14 

23 


X3 
X6 


0.70 ±0.23 
0.84 ±0.28 


1.67 ±0.27 
0.76 ± 0.21 


1.82±»;|f 

i-°i±8:i§ 


27 


X7 


1.88 ±0.37 


1.00 ±0.21 


7-33±g;i 


6-i8±8:H 



ature, Tx ~ 2.2 x 10 7 K, a column density of iV/i = 
1.26 x 10 21 cm -2 and an absorption corrected luminosity of 
L x = (2.34±g;i|) x 10 38 erg s" 1 in the 0.3 - 8.0 keV energy 
band, which are comparable to the values recently obtained 
for the same SNR by Patnaude & Fesen (2003). (After al- 
lowing for differences in assumed distance, they quote a lu- 
minosity of L x = 1.4 x 10 38 erg s" 1 in the 0.5 - 2.1 keV 
energy band, a temperature of T ~ 9.2 x 10 6 K and an ab- 
sorbing column density of Nh = 1.7 X 10 21 cm -2 .) Following 
the same analysis that Blair et al. (1983) performed on this 
SNR, we find firstly that our fitted temperature, is higher 
than the 6 x 10 6 K assumed by Blair et al. (1983, which 



came from the average temperature determined from Ein- 
stein SSS observations of young SNR). This higher tempera- 
ture leads to an age for the SNR of ~ 270 yr and a density of 
120 — 200 cm -3 for the medium into which the SN exploded. 
These results are higher than the age of ~ 120 yr and density 
of ~ 25 cm -3 reported by Blair et al. (1983). Allowing for the 
difference in assumed distances to NGC 4449 would increase 
our age estimate to ~ 380 yr without effecting the density. 
Compared to the non-equilibrium ionization fit for the ion- 
ization timescale, where r = nt — 4.96 x 10 11 cm -3 s, with 
the 90% confidence regions for this parameter giving a range 
of 3.61 x 10 11 <r < 5.00 x 10 13 cm" 3 s, the product of our 
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Figure 6. Left: Hard v. soft hardness ratios for the sources detected in the separate energy bands. The hard hardness ratio is calculated 
from the 2.0 — 8.0 keV and 0.8 — 2.0 keV energy bands and the soft hardness ratio is calculated from the 0.8 — 2.0 keV and 0.3 — 0.8 keV 
energy bands in the usual way as defined in the headings of columns 5 and 6 of Table 3. An indication of the likely origin of the sources 
depending on the location in this plot is also shown by the labels placed next to known sources. Right: Theoretical tracks showing the 
positions different sources would occupy on the hardness ratio plot. The three lines of each type have column densities of 1.4 X 10 20 cm -2 
(Galactic Nh), 1-0 X 10 21 cm -2 (fitted column density for the total emission from NGC 4449) and 3.0 X 10 21 cm" 2 (average fitted 
column density for the fitted sources on the hardness ratio plot) respectively, moving from left to right across the plot. For the three 
thermal models, the temperatures increase from bottom left to top right. Black Body, kT ranges from 0.04 — 1.0 keV; Bremsstrahlung, 
kT ranges from 0.2 — 100 keV; Raymond-Smith, kT ranges from 0.1 — 2.0 keV and the plots shown are for Solar abundance. Decreasing 
the abundance moves the Raymond-Smith tracks to the left. 




37 38 39 40 41 

Log L(X) 



Figure 5. log(TV) - log(Lx) plot for the 24 fitted resolved point 
sources lying within the D25 ellipse. The line shown is the power 
law fit to the high luminosity end of the data. The slope of this 
fit is -0.51 and incompleteness occurs at an absorption corrected 
luminosity of 2.18 X 10 37 erg s~ , which corresponds to a limiting 
absorption corrected flux of 2.12 X 10~ 14 ergs —1 cm -2 . 



calculated values gives 1.0 x 10 12 < r < 1.7 x 10 12 cm" 3 s, 
which gives reasonable agreement for the two methods. 



3.2.4 Source 23 - XRB? 

The position of this source on the hardness plot of Fig. [(J 
places it in the group of XRBs however, its spectrum is best 
fitted by an absorbed non-equilibrium ionization model and 
this source shows no strong evidence for variability over the 
past 10 years, factors which may be more indicative of a 
SNR than XRB. 

3.2.5 Source 27 - XRB 

This source is now the source with the highest individual 
count rate in NGC 4449. Its luminosity has increased by 
nearly an order of magnitude since the ROSAT HRI obser- 
vations of this galaxy in 1994, reported by Vogler & Pietsch 
(1997, see Table 4 for the comparative figures). This is con- 
firmed by the results obtained by fitting it with a 5 keV 
thermal bremsstrahlung model, as used for the ROSAT data,, 
where its luminosity is seen to have increased by a factor of 
at least 3. Its high luminosity and apparent long term vari- 
ability would suggest that it is an high-mass X-ray binary 
(HMXB). The light curve shows no sign of variability within 
our 30 ks observation. 



4 DIFFUSE EMISSION 

A point source and background subtracted spectrum of the 
diffuse emission within the D25 ellipse was extracted and this 
is shown in Fig. |7| The fitted model shown is an absorbed 2 
temperature fit (wabscai (wabs(mekal+mekal))). The com- 
ponent due to Galactic absorption (wabs<3 a ;) has its absorp- 
tion column density fixed at Nh = 1.4 x 10 20 cm -2 and the 
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channel energy (keV) 

Figure 7. Fitted spectrum of the NGC 4449 diffuse emission. 
The fit shown is an absorbed two temperature fit, with a column 
density of Nu = 1.4 X 10 20 cm -2 assumed for Galactic absorp- 
tion and fitted values of N H = (1.29 ± 0.38) X 10 21 cm -2 for the 
absorption within NGC 4449 and temperatures of (0.28 ± 0.01) 
and (0.86 ± 0.04) keV for the two thermal components. The fit- 
ted abundance for the two thermal components was found to be 
(O.32±O.O8)Z 



Table 5. Parameters for the two gas components of the diffuse 
emission. Assumptions: V = 6xl0 65 cm 3 (assumes spherical sym- 
metry); D = 2.93 Mpc; Filling factor, / = 1; n e ~ {El/Vf) 1 / 2 
where EI is the emission integral (normx47rD 2 )/10 — 14 and norm 
is the normalization obtained from the spectral fitting; P ~ 
2n e kT, M ~ n e m p Vf, E Th ~ 3n e kTV and t cool ~ (3fcT)/(An e ) 
where A = Lx /EI. 



Emission 



Soft 



Medium 



kT (keV) 
T CK) 
L x (erg s r ) 
n e (cm -3 ) 
E T h (erg) 
M (Af ) 
P (dyn cm -2 ) 
tcool (yr) 



0.28 ± 0.01 
(3.24 ±0.12) X 10 6 
(6.83±°- 3 |) x 10 38 
0122+ l) ' ()()():1 

U.UIZZICq 0003 

(9.66±° 43 ) X 10 54 
(6.02±°-i|) x 10 6 
(1.07±g-gg) x 10- 11 
(4.66±g-|°) x 10 8 



0.86 ±0.04 

(9.97 ±0.46) x 10 6 

(2.26±°-g) x 10 38 

0059+° ■ 0004 
u - uuoy=c 0.0003 

(1.49±o i 2 ) x 10 55 
(3.01±°- 2 °) x 10 6 
(1.65±g'if) x 10" 11 



(2.02±<]; 3 



10 9 



fit gives N H = (1.29 ± 0.38) x 10 21 cm" 2 for the absorb- 
ing column local to NGC 4449 and values of kT of (0.28 ± 
0.01) keV and (0.86 ± 0.04) keV for the soft and medium 
thermal components respectively. The fitted abundance for 
the two thermal components is (0.32 ± 0.08)Zq. The total 
absorption corrected flux in the 0.3 — 8.0 keV energy band 
for the diffuse emission is (8.86±o;go) X 10 -13 erg s _1 cm -2 , 
which corresponds to an absorption corrected luminosity of 



(9.14±g;§!) x lO 38 erg s . Of this total emission, 75% is 
in the soft component and 25% in the medium component 
and allowance has been made for the flux lost during the 
point source subtraction. The errors shown here are based 
on the 90% confidence levels (1.64er) for the normalization 
values obtained from the fit. Adding a power law compo- 
nent to allow for unresolved point sources within the diffuse 
emission neither improves the fit significantly nor alters the 
fitted parameters for the thermal components. The photon 
index resulting from such a fit was V — 2.53 and the frac- 



tional split of the flux between the 3 components was: soft 
- ~ 70%; medium - ~ 25%; power law - ~ 5%. 

Other parameters of the two gas components have been 
calculated and these are shown in Table 5. The assumptions 
of spherical symmetry and a filling factor of 1 are both likely 
to be over estimates, particularly in the case of the medium 
emission (see below), resulting in the figures being under- 
estimates for n e and P and over-estimates for M, Etk and 
tcool (see Strickland & Stevens 2000 for a discussion of fill- 
ing factors and their likely values in galactic winds). These 
parameters will be discussed further in section 5.2. 

In a further attempt to look at the abundances within 
the thermal components and in particular the ratio of a- 
elements to Fe the spectrum was re-fitted using variable 
abundance mekal models for the two thermal components. 
The abundances of the two thermal components were tied 
and the individual abundances of Mg, Ne, Si and Ca rela- 
tive to Solar values were tied to that of O to form the group 
of a-elements. The fitted abundances for the a-elemcnts 
and Fe were (O.27±o: o |)^ and (O.3O±o:oI)-^0 respectively, 
with a x 2 /d.o.f of 201/117. The a/Fe ratio obtained 
from these values is not significantly different from solar 
((0.91±o'io)-^©) in contrast to the values of Martin et al. 
(2002) for the dwarf starburst NGC 1569, where values of 
2.1 — 3.9Zq are quoted. The difference in these values could 
reflect different contributions made to the hot gas from SNe. 
Type II SNe are the sources of a- elements while Fe is pro- 
duced from Type I SNe. The higher values for the a/Fe 
ratio for NGC 1569 would therefore suggest the presence 
of more Type II SNe than in NGC 4449 and could sug- 
gest that NGC 1569 is experiencing shorter more intense 
bursts of star-formation than NGC 4449 In the latter case, 
as discussed by Delia Ceca et al., 1997, the star-formation 
and energy injection from star-forming regions may be more 
continuous as it is not clear that the present star-formation 
rate is much larger than that in the past on average and 
as such the longer time periods involved will reflect a larger 
contribution from Type I SNe. 

In order to look for variations in temperature within 
the diffuse emission, the region it occupies was divided into 
9 regions and the spectra for each region was extracted and 
fitted. The regions are shown overlaid on the smoothed (us- 
ing a Gaussian with FWHM of 4 pixels ~ 2"), point source 
and background subtracted image of NGC 4449 in Fig. |H| 
and the resulting spectral fits are shown in Table 6. Where 
possible, a 2 temperature fit was used, but for the 4 regions 
with the lowest number of counts (1, 2, 3 and 7) this was not 
possible. Either the fitted temperatures of the 2 components 
were the same, within errors, or the normalization for one of 
the components went to zero. For these 4 regions, the results 
of single temperature fits are given. In addition, for each re- 
gion, the (m — s)/(m + s) hardness ratio has been calculated 
to further highlight variations within the galaxy and these 
are shown in Table 7. To make a comparison of the results 
contained in Tables 6 and 7 clearer, colour-coded maps of 
the temperatures, column densities and hardness ratios are 
shown in Fig. 1101 The two upper panels compare the fitted 
temperatures for the 9 regions. The one on the left shows 
the temperatures of the softer thermal components, where 
two temperatures were fitted and the right panel shows the 
temperatures of the harder components. Where only one 
temperature was fitted, this is shown in both panels. The 
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Figure 8. Low resolution (smoothed using a Gaussian with FWHM of 4 pixels ~ 2") image of the point source and background 
subtracted diffuse X-ray emission overlaid with the 9 regions used to investigate the variation of temperature and absorbing column 
within the diffuse emission. The regions are identified by the overlaid numbers with region 1 in the SW corner and region 9 in the NE 
corner. 



lower left panel compares the column densities (Log(iVf/)) 
fitted to the 9 regions, while the lower right one shows how 
the (m — s)/(m + s) hardness ratio varies. 

These results are not very conclusive but do suggest a 
reduced column density in a region running from north to 
south across the centre of the galaxy and also across the 
northern edge and into the North-East of the galaxy. The 
temperature of the softer thermal components seem to be 
highest in the North of the galaxy and in an arc running 
south across the centre of the galaxy and into the South- 
East. The hardness ratios suggest that generally, the emis- 
sion is slightly higher in the 0.3 — 0.8 keV energy band than 
in the 0.8 — 2.0 keV band. Looking at the soft, medium and 
hard energy band images (0.3 — 0.8 keV, 0.8 — 2.0 keV and 
2.0 — 8.0 keV respectively), a difference can be seen between 
the spatial distribution of the extended diffuse emission in 
these 3 bands, (Fig. [5J . These bands were originally chosen 
to have roughly equal counts in the 3 bands and they high- 
light the relative contributions made by the diffuse emission 
and point sources to the different bands. Surprisingly, the 
0.3 — 0.8 keV band image shows the diffuse emission to be 
extending to the west of the galaxy into regions where the 
spectral fits gave some of the highest column densities. Such 
high absorption should reduce the soft thermal emission and 
suggests that the single temperature fits in these regions are 
not the best models. It is worth noting here that when two 



temperature models were fitted to regions 1, 2 and 3, lower 
column densities were obtained but none of the results were 
well constrained. 



To obtain a clearer picture of the distribution of peaks 
in the X-ray emission, surface brightness slices were taken 
across the galaxy in 4 directions (N to S, NW to SE, W 
to E and SW to NE), centred on the trough of the X-ray 
emission in the centre of the galaxy at a = 12 h 28 m 11.7 s , 
S = +44°06'04.6" with a length of 200" and a width of 20". 
The positions of these slices are overlaid on the Ha image 
of NGC 4449 (kindly supplied by Deidre A. Hunter) shown 
in Fig. 1111 In addition surface brightness profiles were also 
taken from the Ha image using the same slices. The result- 
ing profiles from these slices are shown in Fig. 1121 with those 
from the diffuse X-ray emission, shown by the crosses and 
dashed curves, overlaid on those from the Ha emission (solid 
curves) . Although the X-ray data has been point source sub- 
tracted, the Ha data has not. There is very little correlation 
between the positions of the resolved X-ray point sources 
and the numerous HII regions seen in the Ha data and so 
the slices are being used to look for association of the diffuse 
X-ray emission with the star-forming regions. See Section 
5.1.2 for further discussion. 
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Figure 9. The 3 panels show smoothed images of the soft (0.3 — 0.8 keV), medium (0.8 — 2.0 keV) and hard (2.0 — 8.0 keV) energy 
bands, with the same grey-scale scaling, highlighting the differences in spatial distribution between the bands for the diffuse emission 
and the different spectral characteristics of the point sources. 



5 DISCUSSION 

It is apparent from Fig. that the hardest (> 2 keV) frac- 
tion of the diffuse emission is confined to the central region 
of the galaxy along a ridge running from NNE to the centre. 
This is in the same direction as the orientation of the ma- 
jor axis of the galaxy. The medium (0.8 — 2.0 keV) diffuse 
emission is less extended than the soft (< 0.8 keV) emis- 
sion with the most intense regions of emission in both these 
bands being coincident with the regions where the highest 
density of star clusters is seen (as identified by Gelatt et al. 
2001), suggesting that this emission is associated with the 
increased star-formation occurring in these regions. The less 
intense bulging of the emission seen to the ESE and WNW 
could be indicative of a bipolar outflow along the minor axis 
of the galaxy as would be expected from standard superbub- 
ble models (Weaver et al. 1977). The higher column densities 



seen in Table 6 and Fig llOl for the regions to the SE, SW and 
W of the galaxy would suggest that these areas lie behind 
more absorbing material present in the disk of the galaxy 
and that the eastern and northern edges of the galaxy are 
tilted towards us as shown in the cartoon representation of 
the galaxy's morphology in Fig.^] Allowing for NGC 4449's 
position angle, this would mean the minor axis lies inclined 
to our line-of-sight (tilted both to the N and W) consistent 
with NGC 4449 being a Magellanic Irregular that is inclined 
at 56.2° to our line-of-sight. The lower column density of re- 
gion 3, in particular, would suggest that the X-ray emission 
in that region comes from material lying above the galactic 
disk. 



Chandra Observations of NGC 4449 13 



Table 6. Single or two temperature fits for the 9 regions of the diffuse emission (see text for further details). Columns 1 and 2 contain 
the region numbers and counts for each region respectively. Column 3 lists the column density for the absorbing gas within NGC 4449, 
column 4 the temperature (s) and column 5 the statistic for the fit for each region. Column 6 has the calculated absorption corrected 
luminosities for each of the regions, the errors shown are from the 90% confidence ranges (1.64cr)of the normalizations of the fits or are 
upper limits determined from these values where only one value is given. In all cases, the Galactic absorption column density was fixed 
at N H = 1.4 X 10 20 cm -3 . 



Region Counts N H (cm -2 ) kT (keV) \ 2 / d.o.f L x (ergs -1 ) 

1 197 ±21 2.87 x 10 21 0.19 17.9/14 (2.29± ' 34 ) x 10 38 

2 246 ± 22 3.75 x 10 21 0.21 19.6/15 (3.29± ' 45 ) x 10 38 

3 198 ± 21 4.38 x 10 20 0.45 22.7/14 (3.52± - 70 ) x 10 37 

4 440 ± 26 8.34 x 10 20 0.33,1.03 36.3/20 (1.09±2 |I) x 10 38 

5 1268 ± 39 1.28 x 10 21 0.27, 0.67 55.6/45 (2.98±g :7 |) x 10 38 

6 668 ±30 1.08 x 10 21 0.39,1.13 34.8/29 (1.40±g : ||) x 10 38 

7 45 ± 17 4.69 x 10 21 0.27 4.20/7 (5.44± 47 - 29 ) x 10 37 

8 442 ± 25 1.50 x 10 21 0.23,0.58 19.5/20 (1.33±g;|g) x 10 38 

9 294 ± 23 1.42 x 10 21 0.29,1.93 11.96/17 (8.94± 2; ||) x 10 37 



Table 7. The (m — s)/(m + s) hardness ratio for each of the 9 regions. Column 1 lists the region numbers as shown on Fig. |g] columns 
2 and 3 contain the count rates in the soft and medium bands respectively and column 4 gives the calculated hardness ratios. The softer 
emission is shown by the regions with more negative values. 



Region 


Count rate in 




Count rate in 




(m - s) 




Soft band (counts/s" 




Medium band (counts/s 1 ) 


(m + s) 


1 


(5.33 ±0.57) x 1(T 


■3 


(8.81 ± 0.77) X 10" 


-3 


0.246 ± 0.070 


2 


(7.21 ±0.67) x 10" 


■3 


(2.78 ± 0.49) x 10" 


3 


-0.443 ± 0.091 


3 


(5.39 ±0.59) x 10" 


■3 


(3.26 ± 0.51) x 10" 


■3 


-0.246 ± 0.093 


4 


(11.0 ±0.80) x 10" 


■3 


(7.81 ± 0.70) x 10" 


3 


-0.170 ±0.057 


5 


(31.5 ± 1.20) x 10- 


■3 


(22.4 ± 1.10) x 10" 


■3 


-0.169 ±0.031 


6 


(14.4 ±0.90) x 10- 


■3 


(13.9 ± 0.90) x 10" 


■3 


-0.018 ±0.046 


7 


(1.29 ±0.37) x 10" 


■3 


(0.01 ± 0.38) x 10" 


■3 


-0.985 ± 0.570 


8 


(10.7 ±0.80) x 10" 


■3 


(7.23 ± 0.67) x 10" 


■3 


-0.194 ±0.059 


9 


(5.33 ±0.58) x 10" 


■3 


(6.43 ± 0.64) x 10" 


■3 


0.094 ± 0.074 



5.1 Comparison with Observations at Other 
Wavelengths 

5.1.1 Star Clusters 

The galaxy contains many star clusters spread through-out 
the D25 ellipse and so star-formation is not confined to a sin- 
gle well-defined region within the galaxy. The highest den- 
sity of star clusters identified by Gelatt et al. (2001) is co- 
incident with the peak seen in the diffuse X-ray emission 
at a = 12 h 28 m 11.75" and S = +44°05'36.l". (The largest 
of the contours labelled 3 in the lower panel of Fig. 1131 ) 
The peak in the diffuse emission to the SW of this at 
a = 12 h 28 m 09.65 s and S = ±44°05'11.5" also has a large 
number of star clusters associated with it. It is likely that 
multiple superbubbles are being blown by the combined ac- 
tion of stellar winds and SN explosions from the individual 
stars within these star clusters, leading to the increased dif- 
fuse X-ray emission from hot, shock-heated gas in these re- 
gions. This is supported by the complex morphology seen in 
the patterns of peaks in the X-ray emission shown in Fig. 1121 



5.1.2 Comparison with Ha Emission 

Fig. 1131 shows X-ray contours overlaid on an Ha map of 
NGC 4449 as well as the two separate images. This shows 



that the X-ray and Ha morphologies in the main body of 
the galaxy follow each other closely. The X-ray emission is 
also seen to extend out to the NW into what appears to be 
a hole in the Ha emission, bounded by filaments, labelled 
F in the lower panel. This could be a wind-blown super- 
bubble or the escape of hot gas from a ruptured bubble in 
the form of a super- wind, labelled SB/SW. Ultra-violet ob- 
servations (Hill et al., 1998) have identified regions along the 
north-eastern edge of this hole which show sequential star- 
formation and this activity may be linked to the expansion 
of a super-bubble. As the bubble expands, it sweeps up and 
shock heats the ISM, leading to the formation of areas of 
increased density which can trigger star-formation. In addi- 
tion, the observed filamentary structure seen in Ha could be 
material from a ruptured bubble or entrained ISM material 
caught up in an outflow. A similar, if less pronounced effect 
also seems to be occurring to the SE, labelled SB. The extent 
of both emissions matches well, as can be seen by compar- 
ing the X-ray and Ha surface brightness slices shown for the 
same regions in Fig. 1121 The only exception seems to occur 
to the West of the galaxy where the diffuse X-ray emission 
appears to extend beyond the region of high Ha intensity. 

For the N to S and SW to NE slices, the overall profiles 
of the peaks are similar and they occupy a similar extent. 
The peaks of X-ray emission, in the northern half of these 
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Figure 10. Colour-coded maps of the parameters investigated for variation within the diffuse emission. Top Left: The temperature, in 
keV, of the softer thermal component within the 9 regions. Top Right: The temperature, in keV, of the hotter thermal component, where 
two were fitted. For the four regions with one temperature fits (regions 1, 2, 3 and 7), the soft thermal component temperature is shown. 
Bottom Left: The logged value of the fitted column densities (Log(iV£f )), in cm -2 , for the 9 regions. Bottom Right: The calculated 
(in — s)/(m + s) hardness ratio for the 9 regions 



slices, tend to fall inside the peaks of Ha emission (i.e. they 
are slightly closer to the centre of the galaxy), with the in- 
tensity of the X-ray emission, in the N, being greatest close 
to the Ha peaks. This is the type of behaviour discussed 
by Strickland et al. (2002) in model 5 of the connection be- 
tween X-ray and Ha emission, where the Ha emission is at- 
tributed to a swept-up shell that has cooled, and surrounds 
a hot bubble of SN-ejecta. These potential bubbles in the N 
of the galaxy, labelled B in the lower panel of Fig. 1131 have 
diameters < 1 kpc and the most northern of the two does 
appear to be surrounded by an arc of increased intensity in 
the Ha emission, most clearly visible in the upper left panel 
of Fig. 1131 The increased intensity of the X-ray emission be- 
hind the Ha peaks would be due to conductive evaporation 
of the swept-up material, which increases the density and 



X-ray emissivity of the region where the hot gas and evapo- 
rated material mixes (Weaver et al. 1977). These two slices 
lie closest to the major axis of the galaxy and cut across the 
regions containing the highest densities of star clusters in 
the southern half of the galaxy. The association of the peaks 
is less clear on the south side of both of these slices and this 
may well be a result of the overlapping and merging of many 
wind-blown bubbles from the large number of star clusters 
found there. 

The slices running from NW to SE and W to E lie clos- 
est to the minor axis of the galaxy and cut across the bulging 
of the diffuse X-ray emission. The slices show what appear 
to be cavities in the Ha emission to both the NW and the 
E which are filled with X-ray emission, again suggesting the 
presence of wind-blown super-bubbles, while to the SE and 
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Figure 11. Ha image of NGC 4449 with the positions of the four slices overlaid. Each slice is centred on a = 12 ,l 28 m 11.7 s and 
8 = +44°06'04.7", has a length of 200" and a width of 20". Slice 1 runs from North to South, slice 2 from North-west to South-cast, 
slice 3 from West to East and slice 4 from South-west to North-east. The resulting X-ray and Ha profiles from these four slices are shown 
in Fig.[P2l 



W, the X-ray emission seems to lie outside the Ha emission, 
more suggestive of ruptured super-bubbles, particularly to 
the west as seen by the distribution of X-ray contours, la- 
belled 1, in the lower panel of Fig. 1131 

5.1.3 The Large HI Halo and HI Clouds 

The HI halo of NGC 4449 extends out to a radius of ~ 40 kpc 
(Bajaja et al. 1994) and its outer structure takes the form 
of enormous filaments and clouds (Hunter et al. 1998). The 
VLA observations of Hunter et al. (1998; 1999), show that 
the inner regions consist of a concentration of gas centred 
on the optical galaxy that takes the form of seven large HI 
complexes embedded in a lower density background, with 
the overall outline of the morphology bulging to the ESE 
and WNW as seen in the diffuse X-ray emission. Fig. 12 
of Hunter et al. (1999), shows that the HI in the inner re- 
gions extends beyond the Ha filament seen to the WNW 
in Fig. 1131 and beyond the bulge in the diffuse X-ray emis- 
sion to the ESE. This could be indicative of the swept up, 
compressed ambient interstellar medium lying outside the 
shells of two expanding super-bubbles. The regions of in- 
creased X-ray emission seen in the main body of the galaxy 
as well as lying close to regions of increased Ha emission 
also occupy areas where three of the HI clouds lie. This is 



suggestive of the increased HI density in these regions being 
caused by outflows from star clusters sweeping-up and com- 
pressing the ISM. The overall extent of the inner HI concen- 
tration is ~ 9 kpc and it is embedded in an elliptical region 
of lower concentration with a major axis of ~ 40 kpc from 
which streamers extend. By comparison, the extent of the 
diffuse X-ray emitting region is ~ 2.4 kpc along the major 
axis of the optical galaxy and ~ 1.6 kpc in an approximate 
ESE - WNW direction. The optical extent of the galaxy as 
measured by the D25 ellipse is ~ 4.35 kpc by ~ 3.35 kpc. 
Counter rotation is observed between the HI gas within the 
central region and that outside with the rotation velocity of 
the inner gas being ~ 18 km s _1 and that of the streamers 
seen in the halo ~ 110 km s" 1 . These two velocities and 
the distribution of the HI gas have implications for the po- 
tential of the hot gas, contained within any super-bubbles 
to escape. Modelling the potential of the galaxy as a simple 
spherically-symmetric, truncated, isothermal potential (Bin- 
ney & Tremaine 1987) then the escape velocity at a distance 
r from the centre of the galaxy is given by: 

Vesc(r) = 2 1/2 tv ot [l + ln(r t /r)] 1/2 km s -1 (1) 

where v rot is the maximum rotation velocity of the galaxy 
in km s _1 and rt is the radius at which the potential is 
truncated. If the potential is truncated at the edge of the 
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Figure 12. X-ray surface brightness slices (crosses and dashed curves) overlaid on the Ha surface brightness slices (solid curves). The 
y-axes on the left of each panel gives the counts/degree 2 for the Ho slices while the ones to the right are the values for the X-ray slices. 
Although these axes are on very different scales they do allow a comparison to be made of the spatial distribution of X-ray and Ha peaks 
along each of the slices. All the slices are centred on a = 12 ft 28 m 11.7 s and 6 = +44°06'04.7", located 0.028° along the slice, have a 
length of 200" and a width of 20". This position coincides with the central trough in the diffuse X-ray emission indicated by the small 
central contour labelled 2 at that position in Fig. 1131 



inner HI concentration at r t ~ 4.5 kpc then the escape ve- 
locity for the hot X-ray emitting gas at 1.6 kpc from the 
centre, assuming v rot ~ 18 km s _1 would be ~ 36 km s" 1 . 
This is an upper limit as the value will decrease when the 
extent of the hot gas is deprojected. A similar calculation 
for the potential truncated at the edge of the HI halo at 
~ 40 kpc, with Vrot ~ HO km s _1 gives an upper value for 
Vesc ~ 320 km s . In the absence of radiative cooling, hot 
gas can escape the galaxy's potential if its temperature is 
greater than: 

T esc = 1.5 x 10 5 (uioo) 2 K (2) 

where T esc is the temperature required for the gas to exceed 
the galaxy's escape velocity and uioo is the escape veloc- 
ity in units of 100 km s _1 (Martin 1999). The two escape 
velocities determined above have corresponding escape tem- 
peratures of 1.94 x 10 4 K and 1.54 x 10 6 K respectively. 
From the gas parameters given in Table 5 both the soft and 
medium components of the hot gas have temperatures in ex- 
cess of these values. Eventual escape of the metal-enriched 



SN ejecta contained within the super-bubbles seems possible 
for NGC 4449 but the distribution of the gas in the galaxy's 
extended halo could be as important as gravity in control- 
ling whether the hot gas escapes or not. The presence of 
streamers of HI rather than an homogeneous distribution 
of material can facilitate the escape of hot gas through the 
gaps between the streamers. 



5.1.4 Radio-X-ray Agreement 

The 5 GHz VLA observation (Bignell & Seaquist 1983) 
identifies the known SNR (source 15) in NGC 4449 with 
the most intense peak of the radio emission. The next 3 
most intense peaks on this radio observation also lie close 
to sources in our data (sources 10, 18 and 23) that have all 
been nagged as possible SNRs, although none convincingly. 
In addition, these radio peaks are coincident with the three 
brightest peaks of emission on the Ha map, where intense 
star-formation is occurring and where the presence of SNRs 
in star clusters is to be expected. 
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Figure 13. Top left: Ho grey-scale image of NGC 4449. Top right: Point source and background subtracted diffuse X-ray grey-scale 
image of NGC 4449 on the same spatial scale as the Ha image and with the same smoothing as in Fig. |H| Bottom: X-ray contours 
from the image shown top right overlaid on the Ha image of NGC 4449. The X-ray contours labelled 1 are at a flux density of 
0.65 X 10 -13 erg s -1 cm -2 arcmin -2 . Those labelled 2 are at a flux density of 1.3 X 10 -13 erg s _1 cm -2 arcmin -2 and those labelled 3 
are at 2.6 X 10 -13 erg s _1 cm" 2 arcmin -2 . The background flux density level is ~ 0.25 X 10 -13 erg s" 1 cm -2 arcmin -2 . Also labelled 
are the filaments (F), bubbles (B), super-bubbles (SB) and super-wind (SW) discussed in the text. 
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5.2 Implications of the Calculated Gas 

Parameters on the Fate of NGC 4449 

The presence of what appears to be a developing super-wind 
in this dwarf galaxy can have dramatic effects on its evolu- 
tion. The lower gravitational potential of dwarfs makes the 
possibility of losing newly synthesized metals and swept-up 
ISM material more likely. As such it is worth considering 
the impact that the western super-bubble/super-wind might 
have on NGC 4449. The position of the super-bubble sug- 
gests that its origins may lie in 3 OB associations identified 
by Hill et al. (1994), that lie in the SW region of the main 
body of the galaxy where the highest density of star clusters 
is seen. The ages of these 3 OB associations were found to be, 
5.9, 6.2 and 6.7 Myr while their corresponding masses were 
1.6 x 10 5 , 6.3 x 10 5 and 2.5 x 10 6 M s . The extent that the 
diffuse emission extends from the position of these OB asso- 
ciations is ~ 1.7' or 1.4 kpc for our assumed distance down 
to a flux density level of 0.65x 10 -13 erg s _1 cm -2 arcmin -2 . 
Assuming an average age for the OB associations of 6.3 Myr 
then to have travelled this distance the expansion velocity 
of the bubble has to be ~ 220 km s~ . This will be an un- 
derestimate since the deprojected distance will be greater 
than that assumed and also the presence of a density gra- 
dient in the halo above the disk of the galaxy will lead to 
acceleration of the super-wind. Further, from the calculated 
gas temperature of the soft component of the emission given 
in Table 5, a simple energy conservation analysis would give 
the average speed of a particle to be ~ 280 km s _1 , which is 
in good agreement. At these sorts of velocities, it would take 
~ 200 Myr for the super-wind to reach a distance of 40 kpc 
and so escape the HI halo. This is of the same order as the 
radiative cooling time for the soft component of the diffuse 
emission and again suggests that NGC 4449's huge HI halo 
may allow it to retain its newly synthesized metals and ISM. 
The likelihood that the HI halo will prevent blowout is also 
found when the criterion for blow-out defined by Mac Low 
& McCray (1988) is applied to NGC 4449. This is based on 
the parameter A, defined as the dimensionless rate of kinetic 
energy injection, given by: 

A = 10 4 L m ech,4i^7 p 2 P 4 _3/2 n£ /2 (3) 

where L mec /i,4i is the mechanical energy luminosity in units 
of 10 41 erg s~ , Hkpc is the galaxy scale-height in kpc, P4 is 
the initial pressure of the ISM in units of P/k = 10 4 K cm -3 
and no its initial density. For blow-out to occur, the condi- 
tion A > 100 has to be satisfied. For the whole of the X-ray 
emitting gas, the thermal energy content given in Table 5 is 
2.46 x 10 55 erg. Assuming that this is the result of complete 
thermalization of the kinetic energy of the stellar-winds and 
SN ejecta, occurring over the lifetime of the OB associations, 
then L mec h for the whole galaxy is ~ 1.2 x 10 41 erg s _1 (a 
similar estimate to the 3 x 10 41 erg s _1 obtained by Delia 
Ceca et al, 1997, from the Leitherer & Heckman, 1995, star- 
burst models). The actual thermalization will be less than 
complete, but as argued by Strickland & Stevens (2000), 
will most likely be between 10 and 100 %, meaning that the 
Lmech value above could be greater by up to an order of 
magnitude. Assuming the initial ISM pressure was similar 
to that of our Galaxy then P4 ~ 1 and comparison to the gas 
pressures for the hot gas given in Table 5 show that they are 
an order of magnitude higher which is to be expected for an 



adiabatically expanding over-pressured bubble. As a further 
assumption, the hot gas will be at a lower density than the 
initial ambient density of the ISM and an order of magnitude 
change will be assumed from the average value in Table 5. 
The result of using these figures in the above equation with 
A = 100 means that the scale-height for NGC 4449 has to 
be < 9.5 kpc which is somewhat less than the size of the 
HI halo. In addition this figure will be an overestimate since 
Lmech was estimated from the energy content of the whole 
volume assumed to be occupied by the diffuse emission. The 
volume of the western superbubble is only ~ 10% of the to- 
tal volume and if the ambient ISM density is assumed to be 
closer to that of the hot gas then again the scale-height will 
decrease. If the halo has holes in it, as suggested by the pres- 
ence of streamers rather than an homogeneous distribution, 
then escape is more likely. 

If escape occurs then the amount of mass lost may af- 
fect NGC 4449's ability to continue star- formation. Assum- 
ing the super-bubble has around a tenth of the volume of 
the diffuse emission and as such around a tenth of the mass 
given in Table 5 then over 6.3 Myr, the average mass in- 
jection rate will have been 0.14 Mq yr - , a figure which is 
similar to that predicted for Mrk 33 (Summers et al. 2001), 
another dwarf galaxy that may be developing a super-wind. 
If this rate of mass injection was maintained while the super- 
bubble/super- wind expanded to the edge of the HI halo over 
a time-span of ~ 200 Myr then the total amount of mass 
it would contain would be ~ 2.5 x 10 7 Mq. This value is 
greater than the mass contained within the 3 OB associa- 
tions (~ 3.3 x 10 6 Mq) and the time-scale is also in excess 
of the life-times of the massive stars found there, so unless a 
large amount of mass-loading from ISM material is to occur 
this would be an over-estimate of the mass that the super- 
bubble/super- wind would contain on escape. Compared to 
the total mass of NGC 4449 of ~ 4 x 10 10 Mq this is only 
about 0.1% of the galaxy's total mass. As the time-scale re- 
quired is in excess of the life-times of the massive stars then 
the energy-injection rate would also decrease before escape 
was attained, suggesting that the bubble expansion is likely 
to stall. Hence it seems unlikely that the development of this 
super-bubble/super- wind in NGC 4449 will have a great deal 
of effect on the galaxy's ability to both retain its products 
of star-formation and continue star-formation in the future. 

5.3 Morphology of NGC 4449 

Fig. 1141 shows a cartoon of the possible morphology of 
NGC 4449 and in particular shows how the diffuse X-ray, 
Ha and HI emissions may be associated with each other, at 
the maximum extent of the X-ray and Ha emission, if the 
galaxy is assumed to have a bi-polar outflow in the form 
of two wind-blown super-bubbles along its minor axis. The 
stars represent the approximate positions of the starclusters 
in the two densest populated regions, as identified by Gelatt 
et al. (2001). The structure shown in the X-ray emission in 
the NW super-bubble may be indicative of the bubble rup- 
turing and X-ray emission beginning to extend beyond the 
Ha emission as is suggested in the X-ray and Ha slices of 
Fig. ^| A comparison of the NW-SE and W-E slices which 
run along the northern and southern edges of the bubble 
respectively suggest that the X-ray emission is bounded by 
the Ha emission to the north but not to the South. A similar 
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Figure 14. Cartoon representation of the relationship between the X-ray, Ha and HI morphologies of NGC 4449 and the approximate 
positions of the starclusters in the densest groups (from Gelatt et al. 2001). The extent of the X-ray and Ho emissions were determined 
from a comparison of matched images of the X-ray data and the He* image. The extent of the HI shown is estimated from Fig. 12 of 
Hunter et al. (1998) which shows HI contours overlaid on an image of the same He* data. The structure shown in the X-ray bubble to 
the NW represents the extension of the X-ray emission seen in this direction in Fig. llijl that may actually extend beyond the He* filament 
seen there and as such may be indicative of the escape of hot gas from a ruptured super-bubble in the form of a super-wind. 



picture is seen to the East of the galaxy where the emission 
to the East seems to be bounded whilst that to the SE is 
unbounded. In the regions of intense star-formation to the 
S, there are several peaks of emission in both the X-ray and 
Ha, indicative of a confused morphology where the presence 
of several wind-blown bubbles is to be expected and their 
overlapping both physically and in projection will make indi- 
vidual bubbles impossible to identify. To the N, the presence 
of 2 bubbles with diameters < 1 kpc seems to be indicated 
in the surface brightness slices of Fig. 1121 The HI outline on 
Fig. ll4l represents the extent of the inner HI cloud complexes 
and HI ISM down to a column density of 15 x 10 20 cm" 2 as 
shown on Figs. 10 and 12 of Hunter et al. (1999). Comparing 
these two figures with Fig. ll3l above. although obviously dif- 
ferent in size, all three types of emission do appear to have 
a similar outline to their morphologies. 



6 SUMMARY AND CONCLUSIONS. 

In summary, we have presented an analysis of the X-ray data 
obtained from a 30 ks observation, of the Magellanic Irreg- 
ular, starburst galaxy NGC 4449, by the CHANDRA satel- 
lite. We find X-ray emission from 32 discrete point sources 
in the S3 chip data with 24 of them lying within the opti- 
cal extent of the galaxy, as measured by the extent of the 
D25 ellipse. Some of these can be clearly identified as SNRs, 
others as XRBs and there are also several SSS. We calcu- 
late an age for the previously identified SNR (source 15) of 



270 yr and a density for the medium into which it has ex- 
ploded of being < 200 cm" 3 . The bright SSS (source 14) is a 
candidate for being a white dwarf binary system, while the 
source exhibiting the highest count rate during the obser- 
vation (source 27), shows long term variability and is most 
likely an HMXB. 

The galaxy has a very extended distribution of HI, but 
the inner concentration of HI has an overall general morphol- 
ogy which also follows that of the X-ray emission but on a 
larger scale. The HI complexes within the inner region could 
be the result of outflows from the multiple star clusters in 
the starburst of NGC 4449 compressing and driving the am- 
bient ISM of the galaxy away from the centre, both along the 
minor axis of the galaxy and above regions of increased star 
formation within the disk. The overall extent of the inner 
HI concentration is ~ 9 kpc, that of Ha is ~ 4.35 x 3.35 kpc 
and that of the X-ray emission is ~ 2.5 x 1.6 kpc, with 
this central region being embedded within an huge HI halo 
extending out to ~ 40 kpc. 

The total X-ray luminosity of NGC 4449 in the 0.3 - 
8.0 keV energy band is (2.46±E!i2o) x 10 39 erg s" 1 cm" 2 . Of 
this, ~ 60% is due to resolved point sources and the rest 
is due to diffuse X-ray emission. As suggested by this anal- 
ysis and those of other starburst galaxies (e.g. NGC 253, 
Strickland et al. 2002: NGC 1569, Martin et al. 2002, Della- 
Ceca et al. 1996: NGC 4449 Delia Ceca et al. 1997), the 
diffuse X-ray emission is emerging as a multi-phase en- 
vironment requiring complex spectral models which the 
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spectral resolution of even current X-ray telescopes can- 
not meaningfully constrain. We thus conclude that the dif- 
fuse X-ray emission contains gas at, at least two different 
temperatures and evidence for unresolved point sources. 
The fitted gas temperatures are (0.28 ± 0.01) keV and 
(0.86 ± 0.04) keV for the soft and medium components re- 
spectively and their respective absorption corrected fluxes in 
the 0.3-8.0 keV band are (6.66±g;||) x 10~ 13 erg s" 1 cm" 2 
and (2.20±(i;2|) x 10~ 13 erg s _1 cm -2 corresponding to lu- 
minosities of (6.83±g;||) x 10 38 erg s" 1 and (2.26±g;l£) x 
10 38 erg s _1 . This diffuse emission seems to be more heavily 
absorbed in the SW of the galaxy and shows a higher tem- 
perature in the central and eastern regions consistent with 
NGC 4449 being an inclined Magellanic Irregular galaxy, 
with its minor axis inclined to the N and W of our line-of- 
sight. The morphologies of the X-ray and Ha emissions from 
the galaxy follow each other closely and the X-ray emission 
appears to fill cavities in the Ha emission to the NW and 
E - highly suggestive of wind-blown super-bubbles extend- 
ing along the minor axis of the galaxy. To the SE and W, 
the X-ray emission seems to lie outside the Ha emission and 
this may indicate the developing of a super-wind from the 
rupture of the super-bubbles. The large number of peaks in 
emission seen in both the Ha and diffuse X-ray data indi- 
cates a complex morphology within the main body of the 
galaxy. Regions containing small, (< 1 kpc), bubbles of X- 
ray emission surrounded by shells of Ha emission are seen, 
that are produced by the combined action of stellar-winds 
and SN explosions from massive stars, shock-heating both 
the stellar-ejecta and ISM and sweeping it into shells as the 
ejecta drives its way out of star clusters. 



The hot X-ray emitting gas has a total thermal en- 
ergy content of ~ 2.5 x 10 erg and a total mass of 
~ 9.0 x 10 6 Mq. The super-bubble/super- wind extending 
to the WNW of the galaxy occupies around 10% of the to- 
tal volume occupied by the hot gas (assuming a filling factor 
of 1). The origin of this out- flow appears to be the concentra- 
tion of star clusters observed in the SW region of the galaxy 
disk (Hill et al. 1994). At present, the average energy injec- 
tion rate for the whole galaxy is ~ 1.2 x 10 41 erg s _1 and the 
mass injection rate into the super-bubble is ~ 0.14 Mq yr _1 . 
It seems unlikely that this bubble can escape from the 
huge HI halo at its current estimated expansion speed of 
~ 220 km s _1 , as the time required is in excess of the life- 
times of the massive stars of the star clusters and comparable 
to the radiative cooling time of the soft component of the dif- 
fuse emission, whilst the expansion velocity is less than the 
estimated escape velocity from the HI halo. However, the 
current temperatures of both gas components imply that 
they are capable of escaping the galaxy's gravitational po- 
tential and so venting metal-enriched, hot gas into the IGM. 
The crucial factors are the time for which energy injection 
can be maintained and the actual distribution of HI in the 
halo of NGC 4449. The less homogeneous the distribution, 
the more likely the possibility that some of the hot gas can 
escape. At current mass injection rates, the total mass loss 
would be < 1% of the galaxy's mass and as such would not 
have a catastrophic effect on NGC 4449. 
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